## Loading required package: foreach
##
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
##
## accumulate, when
## Loading required package: iterators
##
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
##
## hour, isoweek, mday, minute, month, quarter, second, wday, week,
## yday, year
## The following objects are masked from 'package:dplyr':
##
## between, first, last
## The following object is masked from 'package:purrr':
##
## transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.040006 -0.003320 0.000519 0.003040 0.120229
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0022478 0.0001405 -15.99 <2e-16 ***
## test_set$FCH4 0.1414673 0.0071422 19.81 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006864 on 2441 degrees of freedom
## Multiple R-squared: 0.1385, Adjusted R-squared: 0.1381
## F-statistic: 392.3 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.155372 -0.002183 0.000236 0.001852 0.097756
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.299e-03 6.117e-05 -21.23 <2e-16 ***
## result$FCH4 4.859e-01 2.776e-03 175.03 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006003 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7582, Adjusted R-squared: 0.7582
## F-statistic: 3.063e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.046687 -0.003331 0.000832 0.003463 0.069299
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0029299 0.0001406 -20.84 <2e-16 ***
## test_set$FCH4 0.1077385 0.0069874 15.42 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006898 on 2441 degrees of freedom
## Multiple R-squared: 0.08875, Adjusted R-squared: 0.08838
## F-statistic: 237.7 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.170860 -0.002059 0.000387 0.002018 0.094157
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.614e-03 6.428e-05 -25.11 <2e-16 ***
## result$FCH4 4.712e-01 2.917e-03 161.52 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006308 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7276, Adjusted R-squared: 0.7275
## F-statistic: 2.609e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.049684 -0.003209 0.000549 0.003025 0.053550
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0023895 0.0001374 -17.39 <2e-16 ***
## test_set$FCH4 0.1334830 0.0067873 19.67 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006744 on 2441 degrees of freedom
## Multiple R-squared: 0.1368, Adjusted R-squared: 0.1364
## F-statistic: 386.8 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.111600 -0.002135 0.000346 0.001931 0.091834
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.492e-03 6.094e-05 -24.48 <2e-16 ***
## result$FCH4 4.741e-01 2.765e-03 171.45 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00598 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7506, Adjusted R-squared: 0.7505
## F-statistic: 2.94e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.057383 -0.003116 0.000761 0.003353 0.091518
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0028336 0.0001358 -20.87 <2e-16 ***
## test_set$FCH4 0.0910863 0.0052743 17.27 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006696 on 2441 degrees of freedom
## Multiple R-squared: 0.1089, Adjusted R-squared: 0.1085
## F-statistic: 298.3 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.166895 -0.002164 0.000562 0.002347 0.079910
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.938e-03 6.801e-05 -28.5 <2e-16 ***
## result$FCH4 4.025e-01 3.086e-03 130.4 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006674 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6352, Adjusted R-squared: 0.6352
## F-statistic: 1.701e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.050944 -0.003289 0.000543 0.003246 0.086515
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0023503 0.0001381 -17.02 <2e-16 ***
## test_set$FCH4 0.1255967 0.0072266 17.38 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006756 on 2441 degrees of freedom
## Multiple R-squared: 0.1101, Adjusted R-squared: 0.1098
## F-statistic: 302.1 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.121750 -0.002120 0.000285 0.001945 0.103527
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.408e-03 6.142e-05 -22.92 <2e-16 ***
## result$FCH4 4.882e-01 2.787e-03 175.15 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006027 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7585, Adjusted R-squared: 0.7584
## F-statistic: 3.068e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.067230 -0.003142 0.000826 0.003442 0.056345
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0027802 0.0001475 -18.85 <2e-16 ***
## test_set$FCH4 0.0932542 0.0066440 14.04 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007261 on 2441 degrees of freedom
## Multiple R-squared: 0.07468, Adjusted R-squared: 0.0743
## F-statistic: 197 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.161268 -0.002152 0.000460 0.002190 0.085538
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.757e-03 6.767e-05 -25.97 <2e-16 ***
## result$FCH4 4.483e-01 3.071e-03 145.97 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006641 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6857, Adjusted R-squared: 0.6856
## F-statistic: 2.131e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.070367 -0.003234 0.000691 0.003299 0.040974
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0024532 0.0001389 -17.66 <2e-16 ***
## test_set$FCH4 0.1285094 0.0061965 20.74 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006815 on 2441 degrees of freedom
## Multiple R-squared: 0.1498, Adjusted R-squared: 0.1495
## F-statistic: 430.1 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.169359 -0.002225 0.000399 0.002048 0.088535
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.487e-03 6.353e-05 -23.41 <2e-16 ***
## result$FCH4 4.518e-01 2.883e-03 156.69 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006235 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7154, Adjusted R-squared: 0.7153
## F-statistic: 2.455e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.041657 -0.003334 0.000564 0.003209 0.067046
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025516 0.0001363 -18.73 <2e-16 ***
## test_set$FCH4 0.1411725 0.0058544 24.11 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006699 on 2441 degrees of freedom
## Multiple R-squared: 0.1924, Adjusted R-squared: 0.1921
## F-statistic: 581.5 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.163063 -0.002234 0.000370 0.002103 0.060342
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.615e-03 6.224e-05 -25.95 <2e-16 ***
## result$FCH4 4.433e-01 2.824e-03 156.95 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006108 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.716, Adjusted R-squared: 0.716
## F-statistic: 2.463e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.058322 -0.003148 0.000593 0.003126 0.045491
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025190 0.0001281 -19.67 <2e-16 ***
## test_set$FCH4 0.1084774 0.0061025 17.78 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006298 on 2441 degrees of freedom
## Multiple R-squared: 0.1146, Adjusted R-squared: 0.1142
## F-statistic: 316 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.148409 -0.002143 0.000404 0.002087 0.078711
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.641e-03 6.185e-05 -26.53 <2e-16 ***
## result$FCH4 4.584e-01 2.807e-03 163.32 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006069 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7319, Adjusted R-squared: 0.7319
## F-statistic: 2.667e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.06683 -0.00382 0.00015 0.00262 0.38233
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0017515 0.0002697 -6.494 1.01e-10 ***
## test_set$FCH4 0.1675527 0.0128210 13.069 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01319 on 2441 degrees of freedom
## Multiple R-squared: 0.06539, Adjusted R-squared: 0.06501
## F-statistic: 170.8 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.19496 -0.00221 0.00009 0.00168 0.37161
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.092e-03 8.816e-05 -12.39 <2e-16 ***
## result$FCH4 5.268e-01 4.001e-03 131.67 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008652 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6396, Adjusted R-squared: 0.6396
## F-statistic: 1.734e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.052202 -0.003838 0.000688 0.003471 0.054760
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0028151 0.0001477 -19.07 <2e-16 ***
## test_set$FCH4 0.0859087 0.0071489 12.02 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007217 on 2441 degrees of freedom
## Multiple R-squared: 0.05586, Adjusted R-squared: 0.05547
## F-statistic: 144.4 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.166321 -0.002195 0.000274 0.001925 0.186303
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.362e-03 7.202e-05 -18.92 <2e-16 ***
## result$FCH4 4.980e-01 3.268e-03 152.39 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007067 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7039, Adjusted R-squared: 0.7039
## F-statistic: 2.322e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.037570 -0.003370 0.000515 0.003148 0.058573
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025736 0.0001416 -18.17 <2e-16 ***
## test_set$FCH4 0.0835142 0.0066356 12.59 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006946 on 2441 degrees of freedom
## Multiple R-squared: 0.06094, Adjusted R-squared: 0.06055
## F-statistic: 158.4 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.191201 -0.002205 0.000319 0.001998 0.099815
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.519e-03 6.668e-05 -22.77 <2e-16 ***
## result$FCH4 4.567e-01 3.026e-03 150.94 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006543 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6999, Adjusted R-squared: 0.6999
## F-statistic: 2.278e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.049401 -0.003553 0.000473 0.003228 0.215192
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0024528 0.0001712 -14.33 <2e-16 ***
## test_set$FCH4 0.1162034 0.0073033 15.91 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008413 on 2441 degrees of freedom
## Multiple R-squared: 0.09397, Adjusted R-squared: 0.0936
## F-statistic: 253.2 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.193797 -0.002251 0.000303 0.002044 0.199596
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.503e-03 7.452e-05 -20.18 <2e-16 ***
## result$FCH4 4.641e-01 3.382e-03 137.24 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007313 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6585, Adjusted R-squared: 0.6584
## F-statistic: 1.884e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.060426 -0.003351 0.000501 0.003240 0.050792
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025748 0.0001328 -19.39 <2e-16 ***
## test_set$FCH4 0.1206494 0.0060113 20.07 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006524 on 2441 degrees of freedom
## Multiple R-squared: 0.1416, Adjusted R-squared: 0.1413
## F-statistic: 402.8 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.112170 -0.002196 0.000342 0.002042 0.076036
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.609e-03 6.172e-05 -26.07 <2e-16 ***
## result$FCH4 4.464e-01 2.801e-03 159.36 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006057 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7222, Adjusted R-squared: 0.7222
## F-statistic: 2.54e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.042367 -0.003556 0.000441 0.003464 0.058800
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0027346 0.0001368 -19.99 <2e-16 ***
## test_set$FCH4 0.0801529 0.0058239 13.76 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00671 on 2441 degrees of freedom
## Multiple R-squared: 0.07201, Adjusted R-squared: 0.07163
## F-statistic: 189.4 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.202201 -0.002323 0.000209 0.002042 0.193362
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.480e-03 7.423e-05 -19.94 <2e-16 ***
## result$FCH4 4.638e-01 3.369e-03 137.69 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007285 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6599, Adjusted R-squared: 0.6599
## F-statistic: 1.896e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.068185 -0.003803 0.000247 0.003104 0.298659
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0023083 0.0002276 -10.14 <2e-16 ***
## test_set$FCH4 0.1388577 0.0100618 13.80 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01118 on 2441 degrees of freedom
## Multiple R-squared: 0.07238, Adjusted R-squared: 0.072
## F-statistic: 190.5 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.163160 -0.002265 0.000230 0.001939 0.288026
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.442e-03 7.978e-05 -18.08 <2e-16 ***
## result$FCH4 4.877e-01 3.620e-03 134.71 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007829 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6501, Adjusted R-squared: 0.65
## F-statistic: 1.815e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.048476 -0.003300 0.000522 0.003111 0.119566
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025336 0.0001343 -18.86 <2e-16 ***
## test_set$FCH4 0.1209069 0.0064598 18.72 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006589 on 2441 degrees of freedom
## Multiple R-squared: 0.1255, Adjusted R-squared: 0.1251
## F-statistic: 350.3 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.173352 -0.002135 0.000328 0.001966 0.094857
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.514e-03 6.244e-05 -24.24 <2e-16 ***
## result$FCH4 4.699e-01 2.834e-03 165.84 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006128 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.7379, Adjusted R-squared: 0.7379
## F-statistic: 2.75e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.043127 -0.003158 0.000817 0.003417 0.036764
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0029084 0.0001279 -22.75 <2e-16 ***
## test_set$FCH4 0.0929055 0.0048287 19.24 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006307 on 2441 degrees of freedom
## Multiple R-squared: 0.1317, Adjusted R-squared: 0.1313
## F-statistic: 370.2 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.193853 -0.002244 0.000617 0.002468 0.076347
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.979e-03 6.651e-05 -29.75 <2e-16 ***
## result$FCH4 3.900e-01 3.018e-03 129.22 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006527 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6309, Adjusted R-squared: 0.6308
## F-statistic: 1.67e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.063779 -0.003421 0.000621 0.003434 0.276040
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0025469 0.0001862 -13.676 <2e-16 ***
## test_set$FCH4 0.0763056 0.0086258 8.846 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009131 on 2441 degrees of freedom
## Multiple R-squared: 0.03106, Adjusted R-squared: 0.03067
## F-statistic: 78.25 on 1 and 2441 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.140648 -0.002291 0.000272 0.001972 0.263238
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.391e-03 7.863e-05 -17.7 <2e-16 ***
## result$FCH4 4.922e-01 3.568e-03 137.9 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007716 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6608, Adjusted R-squared: 0.6607
## F-statistic: 1.903e+04 on 1 and 9769 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.063344 -0.003121 0.000654 0.003023 0.117039
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -0.0024969 0.0001436 -17.39 <2e-16 ***
## test_set$FCH4 0.1065710 0.0061119 17.44 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007044 on 2441 degrees of freedom
## Multiple R-squared: 0.1108, Adjusted R-squared: 0.1104
## F-statistic: 304 on 1 and 2441 DF, p-value: < 2.2e-16




## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.182305 -0.002217 0.000348 0.001957 0.087063
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.511e-03 6.643e-05 -22.75 <2e-16 ***
## result$FCH4 4.379e-01 3.015e-03 145.25 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.006519 on 9769 degrees of freedom
## (19070 observations deleted due to missingness)
## Multiple R-squared: 0.6835, Adjusted R-squared: 0.6835
## F-statistic: 2.11e+04 on 1 and 9769 DF, p-value: < 2.2e-16